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The neutron-neutron fusion process, nn — *> dei>, at very low neutron energies is studied 
in the framework of pionless effective field theory that incorporates dibaryon fields. The 
cross section and electron energy spectrum for this process are calculated up to next- 
to-leading order. We include the radiative corrections of 0(a) calculated for the one- 
body transition amplitude. The precision of our theoretical estimates is found to be 
governed essentially by the accuracy with which the empirical values of the neutron- 
neutron scattering length and effective range are currently known. Also discussed is the 
precision of theoretical estimates of the transition rates of related electroweak processes 
in few-nucleon systems. 



^-mailisandoQtriunif.ca 
2 E-mail: kuboder a@sc . edu 



1 



1. Introduction 

The ultra-high-intensity neutron-beam facilities currently under construction at, e.g., 
the Oak Ridge National Laboratory and the J-PARC are expected to bring great progress 
in high-precision experiments concerning the fundamental properties of the neutron and 
neutron /5-decay. The planned experiments of great importance include the accurate de- 
termination of the neutron electric dipole moment, and high-precision measurements of 
the lifetime and the correlation coefficients for neutron /3-decay. Besides these experi- 
ments that focus on the properties of a single neutron, one might be tempted to consider 
processes that involve the interaction of two free neutrons. An example is the direct ob- 
servation of free neutron-neutron scattering, which would allow the model-independent 
determination of the neutron-neutron scattering length a^ n and effective range r^ n . These 
quantities play an important role in high-precision calculations of low-energy electroweak 
transitions involving the nn channel (see below) . Their accurate values are also important 
in the study of isospin symmetry breaking effects in the strong interactions [1, 2]. Up to 
now, however, because of the lack of free neutron targets, information on Oq™ and r™ has 
been obtained by analyzing processes such as ir~d — > nn^ or nd — ■> nnp; for the current 
status of these experiments, see, e.g., Refs. [3, 4, 5]. Although there has been considerable 
improvement in these analysis, they still contain rather significant theoretical uncertain- 
ties. For instance, the treatments of the exchange currents and higher partial waves in 
the 7i' d — > nri'~f process [6] are open to further examinations; similarly, the treatment of 
the three-body interactions in the nd — > nnp reaction is yet to be settled. Therefore, the 
observation of free neutron-neutron scattering is an interesting alternative, even though 
its extremely low counting rate is certainly a major obstacle. 3 

As another process involving two free neutrons, we consider here the neutron-neutron 
fusion reaction 

n + n — » d + e~ + v e , (1) 

for neutrons of very low energies such as the ultra-cold neutrons and thermal neutrons. 
The experimental observation of this reaction does not seem to belong to the near future 
but, in view of the existing strong thrust for pursuing experiments that take advantage of 
ultra-high-intensity neutron beams, it may not be totally purposeless to make a theoretical 
study of nn-fusion, and this is what we wish to do here. It is worth noting that, for very low 
energy neutrons, the maximum energy E™ ax of the outgoing electron is E™ ax ~ B + 5jy — 
3.52 MeV, where B is the deuteron binding energy and 5jy is the mass difference between 
the neutron and proton, 5 N = m n —m p . This value of E™ ax is significantly larger than the 
maximum energy of electrons from neutron /3-decay, E™p x decay ~ 5^ — 1.29 MeV. Thus 
the nn-fusion electrons with energies larger than Sn are in principle distinguishable from 
the main background of the neutron /5-decay electrons. 

Our study is based on low-energy effective field theory (EFT). For a system of mesons 
and baryons (without or with external probes), EFT provides a systematic perturbative 

3 There exists a project to detect the scattering of two free neutrons using reactor-generated neu- 
trons [7]. 
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expansion scheme in terms of Qj A, where Q denotes a typical momentum scale of the 
process under consideration, and A represents a cutoff scale that characterizes the EFT 
lagrangian [8] . In a higher order effective lagrangian, there appear unknown parameters, 
so-called low energy constants (LECs), which represent the effects of high-energy physics 
that has been integrated out. In many cases, these LECs cannot be determined by the 
symmetries of the theory alone and hence need to be fixed by experiments. Once all the 
relevant LECs are fixed, the EFT lagrangian represents a complete (and hence model- 
independent) lagrangian to a given order of expansion. 

Many studies based on EFT have been made of nuclear processes in few-nucleon sys- 
tems; for reviews see, e.g., Refs. [9, 10, 11] and references therein. These studies have used 
various versions of EFT that differ primarily in the treatment of the following points. (1) 
Choice between the Weinberg counting scheme [12] and the KS W counting scheme [13]; (2) 
Inclusion or exclusion of the pion field as an active degree of freedom [14] ; (3) The presence 
of an abnormally low energy scale reflecting the large scattering lengths in the nucleon- 
nucleon interaction is handled by introducing an expansion of the inverse of the scattering 
amplitude [13], or by introducing the effective dibaryon fields (see, e.g., Ref. [15]). (For 
convenience, we refer to the former as the 1/A expansion method, and the latter as 
dEFT.) For the reasons to be explained below, we use here pionless dEFT. 

The nn-fusion process of our concern here is closely related to the pp-fusion reaction 
[pp — > de + u], and the vd reactions \yd — > ppe(npu), vd — > nne + (npv)], which have been 
studied extensively because of their importance in astrophysics and neutrino physics; see, 
e.g., Refs. [16, 17, 18, 19] for pp-fusion, and Refs. [20, 21, 22, 23, 24] for the vd reactions. 
In an EFT treatment of these reactions, there occurs a LEC that controls the strength of 
contact two-nucleon-axial current coupling. This LEC is referred to as d R in Weinberg- 
scheme calculations [19] and as L\a in KSW-scheme calculations [23]. Park et al. have 
developed a useful hybrid approach called EFT* to determine d R [19] . 4 In EFT*, the 
nuclear transition amplitude for an electroweak process is calculated with the use of the 
relevant transition operator derived from EFT along with the nuclear wave functions ob- 
tained from a high-precision phenomenological nucleon-nucleon potential. 5 An advantage 
of EFT* is that it can be used for light complex nuclei (A=3,4, . . . ) with essentially the 
same accuracy and ease as for the two-nucleon systems. The use of EFT* enabled Park et 
al. to determine d R with good accuracy from the tritium /3-decay rate [19]. 6 This result 
allowed parameter-free calculations of the pp-fusion and f/ep-reaction rates [19], and the 
z/rf-reaction cross sections [24]. It is certainly possible to extend this EFT* calculation to 
the nn-fusion process, but we do not attempt it here. Instead, we employ pionless dEFT 
with motivations to be described below. 

As is well known, the EFT treatment of low-energy nucleon-nucleon scattering is com- 
plicated by the existence of a large length scale associated with the very-weakly bound 
state (deuteron) in the 3 Si- 3 Di channel, or with the near-threshold resonance (virtual 

4 This approach is also called MEEFT (more effective EFT); for a review, see, e.g., Rcf. [11]. 
5 For discussion of some formal aspects of EFT*, see Ref. [25]. 

6 Attempts to determine L\a can be found in, e.g., Refs. [26, 27, 28]; for the possible use of fid capture 
to determine d R or Lia, sec Refs. [29, 30, 31]. 
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state, or singlet deuteron) in the 1 S channel. The 1 /^4-expansion method was proposed 
to cope with this difficulty. This method, however, exhibits rather slow convergence in 
the deuteron channel, for which a typical expansion parameter is (jpd) — 0.4 ~ 1/3, 
where 7 ~ 45.7 MeV and p^ ~ 1.764 fm. Thus, to achieve 1 % or 3 % accuracy, one 
must go up to next-to-next-to-next-to-next-to-leading order (N 4 LO) or next-to-next-to- 
next-to-leading(N 3 LO) order ((1/3) 4 ~ 1 % or (1/3) 3 ~ 3%). These N 4 LO and N 3 LO 
calculations are formidably challenging even in the pionless cases. 7 Another method to 
resolve the difficulty associated with the existence of the large length scale in nucleon- 
nucleon scattering is to introduce explicit S'-channel states, called the "dibaryons" , which 
represent the near-threshold resonance state for the 1 Sq channel, and the deuteron state 
in the 3 S±- 3 Di channel [35, 36]. (As mentioned earlier, we refer to an EFT that includes 
the dibaryon fields as dEFT.) Beane and Savage [15] discussed counting rules in dEFT 
and moreover demonstrated the usefulness of dEFT in describing the electromagnetic 
processes in the two-nucleon systems. Furthermore, a recent study [37] of the electroweak 
observables of the deuteron shows that dEFT allows one to achieve in a very transpar- 
ent and economical way the level of accuracy that, in the 1/^4-expansion method, would 
require the inclusion of terms of very high orders. 

The above consideration motivates us to use dEFT in studying the nn-fusion process. 
We limit ourselves here to the case of very low incident neutron energy, which makes it 
safe to eliminate the pion fields and concentrate on pionless dEFT. We calculate the total 
cross section and electron energy spectrum for nn-fusion up to NLO in dEFT. Although 
the enormous difficulty of observing nn-fusion may render it unwarranted to pursue high 
precision in our calculation, we still wish to aim at a few percent accuracy for the following 
reason. As mentioned, nn-fusion is closely related to pp-fusion and the vd reaction, and 
the latter two processes do require high-precision calculations. Since our present formalism 
can be applied (with essentially no change) to high-precision estimation of the pp-fusion 
and vd cross sections, we may use the nn-fusion case to explain what is involved in those 
high-precision calculations. It turns out (see below) that, at the level of a few percent 
accuracy, we need to treat properly: (1) the LEC, 1±a, which represents the strength of 
a dibaryon-dibaryon-axial-vector (ddA) interaction and which is associated with d R and 
L\a discussed earlier; (2) the radiative corrections; (3) the influence of uncertainties in 
the currently available experimental information on a^ n and rg n . We will show that main 
uncertainties in our calculation of the low-energy nn-fusion cross section come from the 
last item (3). 

2. Effective lagrangian 

For low-energy processes, the weak-interaction Hamiltonian can be taken to be 

H = ^fV M , (2) 

where Gf is the Fermi constant and V ud is the CKM matrix element. / M is the lepton 

7 It has been suggested by Rho and other authors that much better convergence is achieved by adjusting 
the deuteron wave function to fit the asymptotic S'-state normalization constant, Zd = 7Pd/(l — 7Pd) [32, 
33, 34]. 
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current / M = w e 7^(l — 7s)iv, while is the hadronic current, which we calculate here up 
to two-body terms based on the effective lagrangian of dEFT. 

We adopt the standard counting rules of dEFT [15]. Introducing an expansion scale 
Q < A(~ m n ), we count the magnitude of spatial part of the external and loop momenta, 
\p\ and |/|, as Q, and their time components, p° and 1°, as Q 2 . The nucleon and dibaryon 
propagators are of Q~ 2 , and a loop integral carries Q 5 . The scattering lengths and effective 
ranges are counted as Q ~ {7, l/a , l/pd, l/ r o}- The orders of vertices and transition 
amplitudes are easily obtained by counting the numbers of these factors in the lagrangian 
and diagrams, respectively. As discussed below, some vertices acquire factors like r 
and pa after renormalization and thus their orders can differ from what the above naive 
dimensional analysis suggests. 

A pionless dEFT lagrangian may be written as [15, 37] 

C = C N + C s + C t + C st , (3) 

where Cn is a one-nucleon lagrangian, C s is the spin-singlet dibaryon lagrangian including 
coupling to the nucleon, C t is the spin-triplet dibaryon lagrangian including coupling to 
the nucleon; C st describes the weak-interaction transition (due to the axial current) from 
the 1 S'o dibaryon to the 3 Si dibaryon. A pionless one-nucleon lagrangian in the heavy- 
baryon formalism reads 

1 
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where the ellipsis represents terms that do not appear in this calculation, v M is the velocity 
vector satisfying v 2 = 1; we choose v M = (1, 0). is the spin operator 2S' M = (0, <?), while 
Dp = <9 M — \t • V M where V M is the external isovector vector current; A M = — |r • A^, where 
Afj, is the external isovector axial current. g A is the axial- vector coupling constant, and 
m N is the nucleon mass. The terms that involve the dibaryon fields are given by 
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m Ny /r p d _ 

where s a and ti are the dibaryon fields for the 1 Sq and 3 S± channel, respectively. The 
covariant derivative for the dibaryon field is given by = <9 M — iCV^ xt where V^ xt is the 
external vector field. C is the charge operator for the dibaryon field; C = 0, 1,2 for the 
nn, np, pp channel, respectively. a S;t is the sign factor and A s t is the mass difference 
between the dibaryon and two nucleons, m s t = 2m^ + A Sjt . r and p d are the effective 
ranges for the deuteron and 1 S state, respectively. P/ 5 ^ is the projection operator for the 
S = 1 So or 3 Si channel; 



Tr(pf)tpf) = k,, 



(8) 
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Figure 1: Dressed dibaryon propagator (double line with a filled circle) at leading order. 
A single line stands for the nucleon, while a double line represents the bare dibaryon. 




Figure 2: Diagram for the ,S-wave NN scattering amplitude at leading order. The double 
line with a filled circle represents the dressed dibaryon propagator obtained in Fig. 1. 



where Oi (r a ) is the spin (isospin) operator. 

The LECs, y s and y t) represent the dibaryon- NN (dNN) couplings in the spin-singlet 
and spin-triplet states, respectively. These LECs along with A S;i and a Sy t are to be 
determined from the effective ranges in the l So and 3 Si channels. LO diagrams for the 
dressed dibaryon propagators are depicted in Fig. 1. Since an insertion of the two- nucleon 
one-loop diagram does not alter the order of the diagram, the two-nucleon bubbles should 
be summed up to infinite order. Thus the inverse of the dressed dibaryon propagator in 
the center-of-mass (CM) frame reads 



.m N y s 
= %- 



wit \4:7ia s jA S)t 4:7ia Stt E . 

- A — — + — + % V 

4tt m N yj t m N y A t 



(9) 



where we have used dimensional regularization for the loop integral and E is the total 
energy of the two nucleons, E ~ p 2 jm^. The dressed dibaryon propagators are renor- 
malized via the S-wave iViV scattering amplitudes. The amplitudes obtained from the 
diagram in Fig. 2 should satisfy 

4"7T % 

iA s>t = {-iy Stt ) (iD a ,t(p)) (-iy s , t ) = — WAa , Anas t , (10) 

""JV 2 — ^ — U — iy 

where A Sit is related to the S'-wave NN scattering ^-matrix via 

g-l = e"M-l= t f . =-(^U,. mi 

p coto S)t — ip \ 2n J 
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Here S s (S t ) is the phase shift for the 1 S ( 3 Si) channel. Meanwhile, effective range 
expansion reads 

11 1 

pcot5 s = + -r p 2 + ---, pcot5 t = -l + -pdin 2 + P 2 ) + •••, (12) 

CLq z z 

where ao and tq are the scattering length and effective range for the 1 S'o channel and 7 
is the deuteron momentum 7 = y/mjyB (B is the deuteron binding energy) and p d is the 
effective range for the 3 Si channel. Now, the above renormalization condition allows us 
to relate the LECs to the effective-range expansion parameters. For the 1 S , -channel, this 
procedure leads to a s — — 1, 

m N r 1 

d s\jp) = —z — r— 1 — 2 • 13 ) 

m N \ r 2 — + ip - ^r p 2 

For the deuteron channel, one has a t — — 1 and 

2 ft* n / \ m Npd 1 Z d 

yt = — \ — > A(p) = — ^ r- — 1 — 7 2 1 = 7? 1 p ' 14 

m N \ p d 2 7 + ip - ip d (7 2 + P 2 ) E + B 

where Z d is the wave function normalization factor of the deuteron at the pole E = —B, 
and the ellipsis in Eq. (14) denotes corrections that are finite or vanish at E — —B. Thus 
one has [15] 

Z d =-m~. (15) 
1 ~lPd 

This Z d is equal to the asymptotic S'-state normalization constant discussed in Introduc- 
tion. It is to be noted that the order of the LECs y s>t is now of Q 1 ^ 2 , and the deuteron 
state is described by the dressed dibaryon propagator that contains two-nucleon loops as 
well as the bare 3 Si dibaryon. 

The ddA vertex in Eq. (7) contains a LEC Iia, which is associated with the LEC, d R 
or L\Ai appearing in the contact-type two-nucleon- axial- vector vertex. It is not obvious 
how to relate 1±a to d R or L±a, because the dimensions of these LECs are different; Iia is 
a dimensionless quantity, whereas, for instance, L\a is measured in units of fm 3 because 
of two more baryon fields involved in the vertex. However, a relation between Iia and L\a 
is discussed in Refs. [28, 30]. We employ here the assumption proposed in Ref. [37] that 
l\A involves both LO and subleading-order parts. The LO part is fixed from the one-body 
NNA interaction vertex, which is proportional to qa and the factor ( 2^top d ) — 1-024, 
which has been introduced so as to reproduce the result of effective range theory. The 
subleading term Iia represents a two-body interaction and its value is fixed by using the 
ratio of the two- and one-body amplitudes (see below). 

3. Cross section and numerical results 

We calculate the nn-fusion amplitude by adding the contributions from diagrams (a), 
(b) and (c) in Fig. 3. Since the initial two neutrons are in the 1 S state, the dressed 
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dibaryon propagator in diagrams (b) and (c) is limited to be D s {p) (the contribution 
from D t (p) is highly suppressed due to the orthogonality between the initial scattering 
and final bound 3 Si states). Meanwhile, the final deuteron state is described by the 
wavefunction normalization factor \fZ~d- Thus we have the nn-fusion amplitude 
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where e^ d) is the spin polarization vector of the deuteron and e*m is the spatial part of the 
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lepton current Z M = (e^,e*(;)) 



Note that the amplitude in Eq. (16) is proportional to a™ 



and gA- It also depends on 1\a- We remark that the amplitude obtained above is similar to 
that for low-energy np — > d'-f capture calculated in the effective-range expansion approach 
and in an NLO dEFT calculation; the np — > d'y process involves the same partial waves, 
the initial 1 So wave and the final 3 S\ deuteron state. Thus, by changing 



lm N E™ ax in 

the bracket in Eq. (16) to y/m N B, one obtains an expression analogous to the amplitude 
for np — * g?7 (see, e.g., Eq. (39) in Ref. [37]). 

The differential cross section for nn-fusion is now easily obtained. We include here the 
Fermi function, which describes the Coulomb interaction between the out-going electron 
and the deuteron. Furthermore, we take into account the radiative corrections of 0(a) 
calculated for the one-body transition diagrams [38]; here a is the fine structure constant. 
These effects need to be incorporated in order to achieve accuracy better than 1 % in the 
calculated cross section. We then arrive at 8 
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where F(Z, E e ) is the Fermi function defined by F(Z, E e 
27raZ/(3; (3 is the electron velocity (3 = \p e \/E e , and 7„ n = 



x/(l — exp(-x)) with x 
m N E™ ax . Furthermore, 
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where ty is an LEC that appears in calculating radiative corrections. Finally 5^ is the 
radiative correction of 0(a), 9 
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(19) 



withL(x) = / x fln(l-t). 

Eqs. (17) and (18) involve two LECs, ty and l 1A , which need to be fixed. The LEC ey 
can be fixed using the experimental value of the neutron lifetime r and the axial current 
coupling gA- We use the expression for r given in Ref. [38] and employ the experimental 
values, r = 885.7(8) sec and g A = 1.2695(29) quoted in PDG2004 [40]. The LEC, l 1A , 
can in principle be fixed by applying dEFT to the A=3 nuclear systems and using the 
tritium /3-decay rate to constrain l iA , a procedure similar to the one adopted in the 
EFT* calculations [19]. However, a dEFT calculation for the three-nucleon system with 
an external weak current is yet to be done. We therefore make here partial use of the 
results obtained in the EFT* calculations [19]. According to Ref. [19], the cross section 
for charged-current weak-interaction processes in the two-nucleon system receives about 2 
percent corrections from the (two-body) exchange current; see, e.g., Eq. (29) in Ref. [19]. 10 
We may then fix l 1A by imposing the condition that the term involving l 1A should enhance 
the v e d e + nn cross section by 2 % at the initial neutrino energy E v = 20 MeV [41]. 
This requirement leads to 



a 



-^e£ = (2.01 ± 0.40) x 10~ 2 , l 1A = -0.33 ± 0.03 . 



2tt 



(20) 



9 We note that Sa ^ is not exactly the same as the outer radiative correction g{E$,E) given in [39], 
owing to a slightly different renormalization scheme employed calculating 5^ [38]; these two quantities 
are related as <5i^ = g(E ,E) + 5/4. This feature leads to a difference of ~ 1.45 x 1CT 3 in each of 
the subsequent expressions. 

10 Since each scheme of EFT comes with its specific definition of the field variables, the comparison 
of the strengths of the two-body contributions in different schemes involves subtleties. In the present 
exploratory study, however, we simply use the result of EFT* to fix the value of 1\a in dEFT. We will 
give a brief remark on this point later in the text. It is worth mentioning here that the precision of the 
calculated cross section turns out to be controlled by the errors in the currently available experimental 
values of a™ and r^ n . 
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Here we have used Gp = 1.16637(1) x 10~ 5 GeV~ 2 determined from muon decay [40], and 
V ud = 0.9738(4) deduced from the + — > + nuclear /3-decays [42]. The quoted errors in 
ty are dominated by the uncertainties in g A , while the errors in l 1A reflect the ~0.2 % 
error in the calculated v-d cross sections, which in turn are associated with the errors in 
the experimental value of the tritium /3-decay rate. 

We also need to specify the nn scattering length and effective range. Their current 
experimental values are (see, e.g., Ref. [43]) 

a™" = -18.5 ±0.4 [fin], r% n = 2.80 ± 0.11 [fm] . (21) 

Thus there are ~2 % and ~4 % uncertainties in a^ n and r™, respectively. We will discuss 
later the consequences of these uncertainties. We do not include in Eq. (17) the 1/mjv 
corrections or nuclear-dependent (two-body part) 0(a) corrections, the contribution of 
which are about 0.1 %. 

Having specified the LECs and other parameters appearing in our formalism, we are 
now in a position to calculate the cross section and electron energy spectrum for nn- 
fusion. We find it convenient to present our numerical results at a certain specified 
incident neutron energy; the cross sections at other energies can be readily obtained by 
using the 1/v law. So we consider ultra-cold neutrons (UCN). A typical temperature for 
UCN is T UC n ~ 1 m K, and the corresponding average velocity is vucn ~ 5 m/sec. So, for 
the sake of definiteness, we may consider a head-on collision of two neutrons each moving 
with vucn = 5 m/sec; thus in the CM system of these two neutrons v = 2vjjcn ~ 10 
m/sec. The numerical results given below correspond to this kinematics. 

In Fig. 4, we plot the calculated electron energy spectrum, da/dE e , as a function of E e . 
As mentioned, the electrons with E e > 5n = 1-29 MeV can in principle be distinguishable 
from the electrons from neutron /3-decay. Since the amplitude in Eq. (16) is independent 
of E e , the shape of the electron energy spectrum in Fig. 4 is determined mainly by the 
phase factor and the 0(a) corrections coming from the Fermi function and the radiative 
correction, in Eq. (19). 

We also calculate the total cross section a as well as a cut) the latter being the differential 
cross section integrated over E e > 5^; viz. 

p« da r h e da 

a = / dE e —— , a aa = dE e — . (22) 

Jrrie CLJl/ e "'Ojv dlli e 

The results are 

a = (38.6 ± 1.5) x 10~ 40 [cm 2 ], a cut = (30.2 ± 1.2) x 10~ 40 [cm 2 ] . (23) 

Since there is no ^-dependence in the transition amplitude in Eq. (16), the relative errors 
in a and a cut are the same. The ~4% uncertainties in the cross sections in Eq. (23) mainly 
come from the errors in the experimental values of a[J n and rg n , Eq. (21). We note that 
the uncertainties in a^ n and r^™ affect the cross sections to about the same extent; the 
~ 2.2% error in a™ leads to ~ 3.4% uncertainty, and the ~4% error in to ~1.9% 
uncertainty. 
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Figure 4: Spectrum of the electrons from neutron-neutron fusion, nn — > dev. 



4. Discussion and conclusions 

In this paper we studied the nn-fusion process at low energies employing the pionless 
EFT that incorporates the dibaryon fields (pionless dEFT). The electron energy spec- 
trum and the integrated cross sections were calculated up to NLO. We included the 0(a) 
radiative corrections calculated for the one-body transition contributions. Our formalism 
involves the two LECs, e v and 1\a- The former is associated with the inner radiative 
correction in /3-decay, and the latter with the short-range two-nucleon electroweak in- 
teraction. We fixed these LECs with sufficient accuracy for our present purposes in the 
following manner. The LEC, ey, is fixed using the experimental data on neutron /3-decay. 
The LEC, 1±a, is constrained with the use of the results obtained in the EFT* calcula- 
tions in the literature, in which this short-range electroweak effect was determined from 
the tritium /3-decay rate. Once this is done, we can calculate the nn-fusion rate with 
about 4% accuracy. The uncertainties in our theoretical estimates are dominated by the 
existing uncertainties in the measured values of the nn scattering length a™ and effective 
range r™. 

In view of the enormous difficulty of observing nn-fusion, elaborate calculations of its 
cross section are not warranted at the present stage. We have however presented a detailed 
treatment of nn-fusion, because the same formalism can be used for the other related 
processes for which high-precision calculations are certainly needed. The remainder of 
this section is written in the same spirit. 

To fix the value of ey, we have used here the current standard values of Gp , V u d, r, and 
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qa- We note, however, that one of the main purposes of high-precision measurements of 
neutron /3-decay is to deduce the accurate value of V u d avoiding nuclear-model dependence. 
It is therefore important to determine e v through other experiments. 11 Moreover, a recent 
measurement of the neutron lifetime r [46] reported a value that differs from the existing 
world average value by 6.5 standard deviations. A new precise measurement of neutron 
/9-decay currently under planning is essential to clarify this discrepancy and to determine 
the value of V u d- (A change in the rm-fusion cross section due to the new experimental 
value of r is smaller than the uncertainties due to the limited accuracy of Oq™ and r^™.) 

We have fixed the value of 1±a using the result of the EFT* and potential model 
calculations for the tritium /3-decay. It is an important future task to determine 1±a 
within the framework of dEFT itself. For cases that do not involve external currents, 
there exists work in which dEFT is applied to the three- nucleon systems [47]. We need to 
extend this type of work to cases that include external currents. Meanwhile, with the use 
of the value of 1\a deduced in the present work, we can carry out dEFT calculations of 
the pp-fusion process and vd reactions with no free parameters. As mentioned earlier, a 
relation between 1\a and L\a was discussed in Refs. [28, 30]. The relation given in Ref. [28] 
leads to L 1A = -0.26 ± 0.11 fm 3 , while that derived in Ref. [30] gives L 1A = 1.18 ± 0.11 
fm 3 . The indicated errors are consistent with those in the value of L 1A deduced from 
the tritium /5-decay, L\a — 4.2 ±0.1 fm 3 [27], but the central values are smaller. It is 
possible that this discrepancy stems from the different expansion schemes used in these 
calculations. 

We remark that the contribution from the 1\a term probably can play a significant role 
in other processes, such as nd — > nrry and fid — > nnu, which are used for deducing Oq" 
and r™. Recent (pionful) EFT calculations for nd — > nwy [6] and 7c? — > nrm [48] indicate 
that the Kroll-Ruderman (KR) term is important for these reactions. Since the KR term 
is related to the axial-vector couping constant qa through chiral symmetry, and since qa 
is known to be modified by the multi-nucleon effects, it is likely to be important to include 
the short-range two-body effects in analyzing the nd — > nwy and jd — > nnn processes. 
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